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Abstract 

We show in this paper that the 1 /N expansion is a reliable tool to calculate the prop¬ 
erties of a heavy Higgs boson. The 1/N expansion sums up all orders in perturbation 
theory, and therefore avoids the renormalization scheme dependence of the conventional 
perturbative approach. It is explained how effects due to the Landau pole of the Higgs 
sector are isolated and subtracted, and how to perform actual calculations, by comput¬ 
ing the Higgs line shape for the processes ff^H^ ZZ, /'/' at next-to-leading order 
in the 1/A expansion. The results are compared to the perturbative results to show 
the agreement between the perturbative and the nonperturbative approach for Higgs 
masses up to 1 TeV. We conclude that the theoretical predictions for Higgs observables 
are well under control for the entire kinematical region of the LHC. 


1 Introduction 

A main goal of future collider physics is to shed light into the electroweak symmetry break¬ 
ing mechanism. Within the standard model this mechanism is generated by the Higgs 
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sector, leading to gauge invariant mass terms where needed. At the Lagrangean level, this 
is a gauged linear 0(4)-symmetric sigma model, spontaneously broken to 0(3) by a nonva¬ 
nishing vacuum expectation value. The mass of the Higgs boson itself is directly related to 
the corresponding quartic coupling, which implies the nondecoupling property of the Higgs 
boson. For large Higgs mass one expects large quantum corrections which blow up within 
the standard perturbative approach. 

Within the perturbative approach, a lot of work has been done in the recent years. Many 
quantum corrections for relevant observables which contain information on the Higgs boson 
are available |l|. The radiative corrections of enhanced electroweak strength to the decay 
of the Higgs boson were first derived at one-loop in |Q|, and at two-loop in iii- For 
low values of the Higgs mass the theoretical uncertainty is thus small. For heavier masses 
the leading order and next-to-leading order results, start to deviate from the NNLO result, 
inducing theoretical ambiguities. One expects large dependences on the renormalization 
scheme essentially coming from the truncation of the perturbative series which makes theo¬ 
retical predictions unreliable. Whereas the perturbative region will be fully covered by LEP, 
Tevatron and the LHC, it is of special importance to have a good theoretical understanding 
of the heavy Higgs regime which is also in reach of the LHC. Although the existing experi¬ 
mental data points presently in the direction of a Higgs mass in the perturbative region, one 
should not forget that these bounds are logarithmically soft, and that the experimental data 
on sin^^e// still differ considerably in the high precision measurements at LEPl and SLC 

Moreover, certain two-loop quantities are presently known theoretically only as mass 
expansions, and given the high statistics of the experiments, this induces rather substantial 
shifts in the Higgs mass prediction. 

Apart from the phenomenological motivation, it is a major challenge to develop and 
improve calculational methods beyond the standard perturbative approach. A well-known 
nonperturbative approach are Monte Carlo studies on a lattice. The results of lattice 
calculations indicate that for Higgs masses beyond ~ 700 GeV one expects large cutoff effects 

They stem from the fact that the Higgs sector is not asymptotically free. We advocate 
another nonperturbative approach, the 1/A^ expansion. There, the quantum correction 
are ordered by the number of field components |^, |^. The method was applied to the 
Higgs sector to leading order a long time ago |T^ but the large deviation from standard 
perturbative results shaded doubts on its validity. 

Recently we performed a next-to-leading order calculation of the Higgs propagator, 
showing at low coupling a remarkable agreement between the NNLO perturbative result 
and our result |]l^. This proved the 1/N expansion to be a valuable tool for comput¬ 
ing observables beyond the perturbative approach, and for avoiding the truncation of the 
perturbative series and thus eliminating any renormalization scale dependence. 

In this paper we extend the results to the case of three-point functions, which results 
in a quantitative, nonperturbative understanding of the physical processes // ^ H —> 
f f, Z^Zl. It was tried to make the computation as transparent as possible to convince 
the reader that the method is conceptually well settled and practicable, though some amount 
of numerical work is unavoidable, which on the other hand is common to the calculation of 
higher-order finite momentum Green functions in massive theories. 

After some preliminaries, where we remind the LO results and the issue of the tachyon 
of the 1/N approach, we show how one obtains a well defined, tachyon free representation 
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of Green functions. An ultraviolet subtraction scheme is defined in terms of an order 1/N 
counterterm Lagrangean. Then, the finite and tachyon free Green functions are expressed 
in form of multiloop graphs made out of dressed propagators. These can be reduced to two- 
dimensional integrals over dressed propagators and form factors. Then we describe how to 
perform the remaining integrations numerically, to obtain the final results for the two- and 
three-point functions. We show and discuss our nonperturbative result for the Higgs line 
shapes in two physical scattering processes, and compare to the perturbative result. After 
this discussion, we summarize the conclusions of the paper. 

2 Preliminaries 

In this section we review the leading order result for fixing our notations, and discuss 
the structure of the counterterms which absorb the ultraviolet divergencies at leading and 
next-to-leading order. 

Because we are interested in the effects of a heavy Higgs boson, we neglect all the gauge 
couplings and also the higher order corrections in the Yukawa couplings. Thus, the starting 
point for the 1/N expansion in a scalar theory is given by an 0(A^)-symmetric (j)^ theory: 

^ = ( 1 ) 

Here (po is a Ai-component vector and no, Aq are bare coupling constants indicated by the 
subscript zero. Because we want to focus on the spontaneously broken case, we rewrite it 
as 


C = ^ ( 2 ) 

Nvq = —QN^q/Xq is the bare vacuum expectation value of (j)^. The symmetry of the 
ground state is broken down to 0{N — 1), which amounts to the presence of A^ — 1 massless 
Goldstone bosons (GBs), 7ri=i^...^Ar_i, and a Higgs particle h, of mass = AqUq/S. In the 
case Ai = 4, the GBs are related to the longitudinal degrees of freedom of the vector bosons 
Wl, Zl by the equivalence theorem. 

To overcome the combinatorial difficulties of multi-loop calculations, we use a well- 
known auxiliary field formalism |^]. We add to the Lagrangean the following expression, 
which becomes identically zero when one uses the equations of motion: 


Tv = 


3N 




Xo - -^{^0 - 


(3) 


is an arbitrary Lagrange multiplier. One additionally has the freedom to replace xo by 
axo + /?) where a, (3 are some constants, without changing the equations of motion. These 
constants correspond to a x vacuum expectation value and a wave function renormalization. 
We introduce renormalized quantities in the following way: 
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</’o — 




(4) 


From the discussion above it is clear that only Z/j, Z\ and Z^ are renormalization 
constants in the usual sense. The other can be chosen arbitrarily as long as the choice is 
consistent. Their introduction will turn out to be useful for the renormalization of functions 
which contain auxiliary particles as external legs. The Lagrangean is now of the form: 


r , Ar 2 C2 2 

^ = —dy(^d(^2cxNvx2— — 2 “^ 

2 2 

- iV..= ) - - Nv--) 


with 


8 Nv^ ^ 


(f) = (yf, C6 (T + a//Vu) 


(5) 


Note that a mixing between y; and a takes place. We see that the expression contains 
— up to a trivial surface term Dy — all possible field combinations of dimension four 
which are built out of an 0{N) vector and a dimension two singlet field allowed by Lorentz 
invariance and the 0{N) symmetry. The vacuum expectation value of the Higgs reduces 
the symmetry as usual to 0{N — 1). The constants Ci=o,...,6 are expressible in terms of the 
above renormalization constants. 


CO = Z,^ Cl = 

C3 — Z^Z^ C4 — Z^ZfjV^ 

C6 ~ \J^h/ 


C2 


Zx 


C5 = (1 — Z«;)Za 


( 6 ) 


These constants define a redundant counterterm structure of the theory, and are fixed 
by imposing conditions on one-particle irreducible Green functions. Whereas some are 
necessary to absorb the divergencies occurring in physical processes, other can be adjusted 
to make subtractions on Green functions containing auxiliary fields as external legs. Clearly 
all these subtractions have to be consistent with each other. In order to determine the 
renormalization constants, the quantum corrections to the tree level Green functions have 
to be calculated up to a given order, and have to be related to the counterterms by a 
renormalization scheme. 

In perturbation theory the occurring renormalization scheme dependence is proportional 
to some power of the coupling constant, and namely one order higher than the loop order 
in which the process is calculated. Thus it is enhanced in the case of strong coupling. It is 
a crucial point that this is not the case in the 1/N expansion. Once the relations between 
a set of observables and the Green functions are fixed, all further predictions are exact, up 
to unevaluated powers of 1/N. A change of the ultraviolet subtraction point fj? (see below) 
is compensated by a change in A. Thus it does not effect any observable prediction, being 
merely a reparameterization of the theory. 
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Figure 1: Leading order tadpole and self-energy graphs. 

Before describing the computation of the Higgs propagator in next-to-leading order in 
1/A^, we would like to repeat the leading order results to set the stage for the next step. To 
order (1/iV)® the quantum corrections are simply the sum of all graphs containing arbitrary 
numbers of external y legs attached to a GB loop. The sum can be carried out and a 
simple closed expression for the resulting effective potential can be found |^. The only 
divergent graphs are a y tadpole and a momentum dependent x self-energy graph shown 
in figure 1. Both divergencies can be absorbed into and Zx, or in other words, into the 
coupling constants of the theory. No other renormalization constant is involved. Whereas 
the observable vev is defined by the ground state of the theory — resulting in an observable 
vector boson mass if one couples the model to a gauge theory — the value of A can be fixed 
by a measurement of the Higgs mass. In terms of the introduced renormalization constants, 
one has; 

Zl = l-2T^, 1/Za = 1 - AH0( V)/3, Z^ = Zh = Z^ = Z^ = l, = 0 . (7) 

The integrals which correspond to the leading order corrections are trivial. Absorbing 
combinatorial factors in the expressions, one finds in dimensional regularization: 

.ri/, 1 , const. 

-4°(s) = -^^log(-s/A^) H-^-h const.’ 

= 0 ( 8 ) 

Defining an energy dependent Higgs mass by m?{s) = u^A(s)/3, where 1/A(s) = 1/A -|- 
(q;o(s) — ao{p‘^))/3, one gets an expression for the matrix of two-point functions. Its inverse 
defines the dressed propagators to leading order, /r is the renormalization scale. The 
resulting propagators are found to be: 


Dhh 

1 

s — m?‘{s) 

Dhx 

i m‘^{s) 

y/Nv S - ■W?‘{S 

^xx 

i sm?{s) 
Nv"^ s — m?{s) 

D-lTlTTk 

= -6ik 
s 


The GB propagator is not influenced at leading order. Glearly the same results can be 
obtained without the auxiliary field. The dressed x propagator corresponds to a geometric 
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series of loop graphs E3- Note that the propagators contain an imaginary part for physical, 
timelike momenta, which describes Higgs decay into Goldstone bosons. To keep the pole 
on the correct Riemann sheet one has to use the ie prescription in the logarithm of the 
s-dependent mass: 


1 

m?{s) 



1 

327r2u2 


log( 



( 10 ) 


For spacelike momenta, a Landau pole occurs at the scale: 


Al = ^ exp 



( 11 ) 


Since m?{s) is negative beyond the Landau scale A^, all these propagators have an extra 
pole for a spacelike momentum s = which is defined by the equation: 

T ?7i^(—Ay) = 0, Ay > 0 . (12) 


There is a tendency to discard non asymptotically free theories as fundamental theories. 
The occurrence of the tachyon in the sigma model is widely believed to be an indication of 
an ill-defined quantum field theory. However, to or knowledge there is so far no rigorous 
proof thereof. We will show in the next section that the tachyon is technically an artifact 
of perturbation theory and can by no means be considered an inconsistency of the model. 
Apart from these issues, it is natural to ask for higher order effects. How to go beyond the 
leading order in 1 /N is directly related to the tachyon problem, because on has to integrate 
over graphs built out of tachyon contaminated propagators. A method to treat this is given 
below. Such a procedure is touching the aspects of the theory which cannot be investi¬ 
gated by Feynman diagrams, even summed up to all orders. In this sense any tachyonic 
regularization is modeling either effects from a more fundamental theory containing the 
Higgs sector as a low energy limit, or a certain type of nonperturbative effects within the 
Higgs sector itself. Our nonperturbative treatment of the sigma model is independent of an 
interpretation either as an effective theory or as a fundamental theory, and we will make 
no assumption in this direction in this paper. 


3 Calculation 

In this section we explain the necessary steps to calculate the NLO corrections to the pro¬ 
cesses // —> iL —> ZlZl^WlWl, f f in the 1/N expansion. In the case of longitudinal 
vector boson decay modes, we calculate the processes by using the equivalence theorem. 
Therefore we repace the longitudinal vector bosons by the corresponding would-be Gold- 
stone bosons z and w. To this end one has to know the two- and three-point functions 
to that order. To get a physically meaningful quantity one has to define a perturbative 
renormalization procedure to deal with the usual divergencies. On the other hand one has 
to define a scheme to treat the tachyon which is a perturbative artifact of the leading or¬ 
der calculation. We have seen that the infinite series of Feynman diagrams shows up in 
the propagators as a momentum dependent mass producing a tachyon pole. The dressed 
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propagators in eq. ^ are, together with the occurring vertices and the counterterms dehned 
in eqn. |^, the building blocks of the diagrams which describe the next-to-leading order 
effects. To leading order, only GB propagators appear as internal lines. It is crucial for the 
simplicity of the calculation that, with the help of the auxiliary field Xi the next-to-leading 
order contributions of the 1/N expansion to the Higgs propagator can be depicted by a 
small number of graphs. 

In the following, a tachyonic regularization will be defined. Then the diagrammatical 
part of the calculation will be given together with the renormalization procedure. The 
subtractions will be explained in detail for the relevant multi-loop graphs, and the finite, 
tachyon-free quantum corrections to the two- and three-point functions will be derived. 
They will be represented as two-fold integrals. Finally the remaining two dimensional 
integration is performed numerically, which gives quantitatively the propagator and vertex 
corrections. 


3.1 Tachyon subtraction 


Before starting the calculation, we introduce a tachyon-free representation for the occurring 
propagators |11]. Assuming that its occurrence is not due the failure of the theory under 
consideration, but of the intermediary expansion technique used |n], it is reasonable to 
circumvent the ill-defined part by an adequate subtraction of the tachyonic pole. If the 
theory is only describing a low energy effective theory, some physics at a higher scale will 
cure the problem. It is clear that a priori the tachyon regularization is not unique, and one 
has to pay the price in the form of an additional ambiguity. We modify the propagators 
in eq. ^ by the following subtraction of the tachyonic pole, which is the minimal way of 
eliminating the tachyonic pole: 


1 


1 


1 


s — m?[s) 
m^(s) 
s — m?{s) 


s — m?{s) s + 
m?{s) Aj, 

s — m?{s) s -|- A|, ^ 


(13) 


K = 327r^u^/(327r^u^ -|- Aj.) is just the residuum of the tachyonic pole. To get a quantitative 
feeling for the size of the subtracted terms we give some values of the appearing constants 
in table ||. It is clear from eq. n that only a pair of parameters m/i,// defines a physical 
situation. We related p to the position of the Higgs pole in the complex plane, Mp^ie — 
iTpoiel‘2', by the relation + rp^jg/4 ||T^. Its position is also shown in the 

table. Near the Higgs resonance, the residuum of the tachyonic pole is a measure for the 
quantitative influence of the tachyonic regularization. Any amplitude will be modified there 
by terms proportional to k. The perturbative expansion of these terms vanishes identically. 

Another strategy to treat the tachyon issue would be the introduction of a momentum 


cutoff IQ . Such a cutoff has to be under the tachyon scale, and destroys the gauge invariance 
in the gauged model. On the contrary, our tachyonic regularization preserves the symmetry 
of the model in all orders of perturbation theory. 


7 











mh 

Mpole 

^ Pole 

Al 

Ap 

K 

0.20 

0.20 

5.3x10"^ 

1.7x10^^ 

1.7x10^^ 

1.6x10"^'^ 

0.40 

0.40 

0.04 

1.2x10^ 

1.2x10^ 

0.1x10““ 

0.60 

0.59 

0.15 

0.5x10^ 

0.5x10^ 

0.6x10“^ 

0.80 

0.74 

0.34 

3.2x10 

3.2x10 

4.7x10“^ 

1.00 

0.83 

0.58 

9.64 

9.88 

4.7x10“^ 

1.20 

0.87 

0.82 

5.03 

5.45 

0.14 

1.40 

0.86 

1.00 

3.36 

3.93 

0.24 

1.60 

0.84 

1.11 

2.57 

3.23 

0.31 

1.80 

0.82 

1.19 

2.12 

2.85 

0.37 

2.00 

0.81 

1.24 

1.85 

2.62 

0.41 


Table 1: Information on the physical and the tachyonic pole at leading order. The param¬ 
eter is related to the physical pole in the way explained in the text. The dimensionful 
quantities are in TeV. 


^TTTT 
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MwmNvx- 
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AAA/»- 

^XX — 
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Figure 2: Feynman rules. The dot indicates the dressed, tachyon free propagators. The 
counterterms are defined by eg. 14 - 


3.2 Diagrammatics 

The Feynman rules for the next-to-leading order calculation contain the full LO propagators 
given above in eq. and the interaction terms defined in eqs. S. The propagators are un¬ 
derstood as tachyonicaly regulated by the prescription of eqs. Thus the diagrammatical 
part of the calculation is settled. The Feynman rules are given in figure 

The dot indicates the dressed propagators stemming from the leading order calculation. 
We now fix the counterterms. Considering the values of the leading order renormalization 
constants, we write = 1 — k/N, Z,p = 1 z^/N, Zh = 1 + Zh/N, Z^ = 1 2zylN, 

= v^/N, XjZx = 1 — zx/N , Zj^ = 1 + z^jN . The LO counterterm is already absorbed in 
the definition of the dressed propagators. The NLO counterterm Lagrangean is expressible 
in terms of these new variables: 
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Figure 3: Relevant two- and three-point functions to be calculated. 


C, 


counterterm — 7 - 


Nv'^i^ + Zv)x + (2z^ - K- zx)^ 

2 ^ 2 

Zy — K Ar„.2\ Ar„,2\ ^^2 Ar„.2\2 


= (^, 





+ \/ivu) 


((/)^ - 


(14) 


As we will see below, the NLO corrections are expressible in form of multi-loop diagrams 
made out of the dressed propagators. The counterterms in eq. 14 will allow to make 


subtractions of subgraphs inside these diagrams which is the motivation for our redundant 
renormalization formalism. 

The relevant two- and three-point functions are shown in figure Their actual form in 
terms of graphs will be shown below. 

The counterterms are fixed by some renormalization conditions. Two conditions are 
already dictated by the cancellation of tadpole contributions. Further conditions may be 
imposed on one-particle irreducible two- and three-point functions. We consider the follow¬ 
ing conditions: 


r7r7r(0) ~ ^hhi^) ~ ~ ~ ~ ~ r/i7r^(0, 0) — 0 (15) 

These conditions define a renormalization scheme. Otherwise, all possible counterterms 
are fixed by the symmetry of the counterterm Lagrangean of eq. E- Whenever possible 
we set the renormalization scale to zero, because the multiloop counterterms turn out to 
be analytically simpler, as will become clear later. When the subtraction point // cannot 
be set to zero because of infrared problems, we combined the multiloop graph with the 
corresponding counterterms at an Euclidian p and such that the sum is free of the infrared 
singularity, and afterwards we took the limit ^ ^ 0. 

The list of needed graphs is given in figure Q. One finds these graphs by usual power 
counting in 1/A^. A detailed discussion of their computation will follow in the next subsec¬ 
tion where the analytical part of the calculation is done. The renormalization constants can 
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Figure 4: List of occurring graphs to order 1/N 
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be determined through the scheme of eq. El in terms of vacuum graphs of the type given 
in figure By using 


i(/i2) = _i,(_^2)/2 + i2(0)/2-i3(0)-i4(0)/2-i5(0) + i6(0)/2 

— As{—ii‘^) + Ag{—^‘^)/4 
C(^2) = ^^(0)_C2(0) + C3(V) 

one finds 

z^ = -v^D'{0) Zh = -v‘^B[{0) Zy = -f^-^^ 

^x = S^2(0) + C(0) «=^^2(0) v^ = v‘^n ' (16) 

z^ = ^A{^^^) + ^^BM + 2C{^A^) 

h 

Note that the occurring expressions can be expanded in A. This provides the connection 
to the standard loop expansion, and allows to check the consistency of the renormalization 
procedure to any order diagramatically. 


3.3 Analytical framework 

In order to keep our expressions compact, we introduce the following notations: 


9{s) 



fis) 


m?{s) 

S — 771,2 (s) 


(17) 


We further do not write down explicitly the tachyonic regularization of the propagators, 
but we understand that the tachyonic pole has to be subtracted according to the discussion 
of the previous section before numerical integration is performed. 

To obtain compact formulae for the occurring 1/A graphs, we use in the following the 
massless scalar one-loop three- and four-point functions as building blocks. Their explicit 
expressions []l^ |l^ are listed in the appendix. With the help of these functions, the two- 
and three-loop topologies can be cast into the form of two-fold integrals. The 1/N graphs 
which we need correspond to the following list of dimensionless integrals, where r = (q + k) 
and we use the notation dQ = k and /r stand for a timelike and spacelike four- 

vector respectively, and is a lightlike four-vector corresponding to on-shell GBs. To work 
with dimensionless quantities, we divide the 1/A graphs by appropriate factors of v. We 
distinguish ultraviolet subtracted graphs from unsubtracted ones by a tilde. 


3.3.1 One-loop topologies 


fh 


^ J dQf{q^) 

- (^T^) ^ 


(18) 

( 19 ) 
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^2 

Bi 

B2 

Cl 

D 

E 


Fi 


F2 


j dQg{q^) gir"^) - g{q^ 

^ J dQq^f{q^) gir"^) - g{q^) - 
^ J dQf{q^){f{r^) - f{q^)) 

X [ dQf{q‘^)\g{r'^) - g{q^ 


d 


dQq^f{q^ 


-1) dQ 


{q + iy 


dQf{q^ 


1_1_,2A| 1 

r2 g2 ^J.2\k^=0^2 

f{q^)f{{q-kf)-f\q^) 

q^ 1 ^ 


(q — iY{q — I + k)"^ (?2 


dQ 


;f{q^)f{{q-kf)-f\q^) 


{q + iy 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 


All necessary counterterms for the subtractions appear in the counter Lagrangean of 
eq. ( 0 ), as one can show by using the relations Th = D{0) = -Bi(O) — - 62 ( 0 ), Fi(0) +-F 2 ( 0 ) = 
( 71 ( 0 ). For convenience we define Bi{s) = Bi(s) — -Bi(O) — sD'{0) because we will not 
explicitly need the wave function renormalization of the Higgs boson which is related to B'. 
Note that T/j = -D(O) is a necessary consistency condition to keep the GBs massless after 
summing the perturbative series. This holds independently of the tachyon. One can check 
that to render Bi and D finite, one subtraction is already sufficient. This means that one 
gets an ultraviolet convergent wave function renormalization constant, in contrast to usual 
perturbation theory. Also note that graphs like B 2 , with two momentum dependent masses 
in the numerator, have an extra suppression factor of l/log‘^{q^) in the ultraviolet, which 
leads to a convergent integral^. The summation of all loop graphs thus improves the UV 
behaviour. This ultraviolet screening is a remarkable feature of the nonperturbative \/N 
expansion. 


3.3.2 Higher-loop topologies 


For the higher-loop topologies one has to subtract the subgraphs first. The necessary 
subtractions are most transparent in a graphical representation. We have collected them in 
figure 

The reader can check that the counterterms defined in eq. 


14 correspond to sums of 


boxed graphs. For the three-point counterterms, this is shown in figure g. 

Note that the sum of the 1/N graphs and counterterms for a given topology gives an 
infrared finite result. Thus, to derive the following formulae, it is necessary to introduce an 
infrared regulator for the individual terms which cancels in the final sum. As stated above. 


®The integral converges as 1/ log Acutof f 
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Figure 5: Subtraction scheme of the higher loop topologies. In the countergraph of D, the 
dash stands for differentiation with respect to the external momentum and the box means 
that the contained subgraphs have to be evaluated at vanishing external momenta. 



Figure 6; Graphical representation of the F^^jr counterterm. 
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we are using analytical expressions for the massless three- and four- point functions, Cq, 
-Dq) which are given in an appendix for completeness. 


^3 


A 4 

^5 

^6 


C 2 


F 3 


^ J dQdP 




(p 


,2^2 


1 

11 

_{p + ky 

q2 


1 


\ y 


(p + q)'^ p^ {p^Y 


V 


kf'iQ— 


1 


dQ 


g2(47r)2 

f{q^ 


1 - 


{q^fDo{k,q,-k) - 


2{q-hf q^ 

q2k2 ^Vg2 


1 + logd?) 


Stt^ 


dQ 


dQ 


f{r^)f{q^)Co{k,q) - f{q^) 


2 -|- log(p 

(47r)2(72 


r\^f{r^)f{q^)Ci{k,q)-f{q^) 


2/„2x 


2 + 


i°g(g) V 


(47r)2 I 


^ / dQ 


-I dQ 


q^f{r^)f{q^)Co{k, q) - f{q^)^ ^ 


(47r) 


_/2^^2^2 + log(|2 


f{r‘^)fiq‘^)Co{k,q) 


{q-if 


(47r) 


p2gi2 


(29) 


We will need in the following sums of these expressions taking into account combinatorial 
factors and signs. 

A{s) = --yli(s)-I--^ 2 ( 5 ) - 213 ( 5 ) -- 214 ( 5 ) - yl 5 (s)-I--y46(s) 

B{s) = B 4 {s) + B 2 {s) 

C{s) = Ci(5)-C2(5) 

F{s) = Fi(5) + F2(5)-F3(5) (30) 


The functions 2 I, B, C, D, E, F are integrals over Minkowski space, parameterized by qi, = 
{qQ,pn). After performing the angular integration, the number of integrations can be re¬ 
duced to two. Choosing the timelike momentum as ky = (-^/s, (l), it is obvious that the 
integrals which belong to two-point functions are isotropic. For the three-point function 
with on-shell GBs, one has a second timelike four-vector ly = (1, 0,0, l)y^/2. Explicitly, 
one finds for the angular integration; 


dkl 


1 


1 


dQ. 


dvr {q — ly 

1 


dvr {q — iy{q — l + kY 


2py/s 

1 

Aqops 


log 


log 


q^ + \/g(go + p) 
q^ + \Qs{qo - p) 

{q^y - s{qo - pY 
{q^y - s{qo + pY 


(31) 
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Figure 7: Integration contours for the numerieal integration of the (a) two- and (b) three- 
point funetions. e = y/s. 


Thus we end up with two-dimensional integral representations for the relevant graphs con¬ 
tributing to the processes we look for. In principle the qq integration is also manageable 
by complex rotation, but one has to take great care because of the complicated singularity 
and cut structure of the integrands. We proceed numerically in the following. 

3.4 Numerical calculation of the diagrams 

The numerical problems encountered in the calculation of the needed diagrams are twofold. 
First of all, one faces convergence problems due to the physical particle poles. Although 
the massive propagators have a finite width in our case, one still has to use a Wick rotation 
to get a fast numerical convergence. To do the two dimensional numerical integration: 


OO OO 

J dqo j dpp^f{qo,p) 

—OO 0 

one chooses the contours in the complex p plane to avoid the cuts of the logarithms and to 
stay away from the poles. They are shown in figure |^. For the two-point functions, a simple 
rotation p —> /9exp(—i(/)) leads to a smooth integrand (a), for the three point function one 
has to avoid also the cut of the logarithms present in eq. ^ (b). No singularity and cut is 
crossed for any f G (0,7r/2). The independence of the numerical result of f serves as a test 
of the stability of the integration. 

The procedure gives fast and accurate results as long as the integrand is well behaved 
at infinity. This is the second problem of the numerical part of the calculation. Usually in 
loop calculations one is able to separate the UV divergent part by evaluating the Feynman 
graphs in 4 -|- e dimensions. The regularization leads to manifestly finite integrands. If, as 
in our calculation, the regularization has to be done by subtracting divergent counterterms 
from divergent graphs, the numerical situation is more delicate because of cancellations. To 
ensure numerical ultraviolet stability, it is necessary to arrange the integrand in such a way 
that as many cancellations as possible are manifest. For the integrations we used adaptive, 
deterministic Fortran routines. 
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4 The results 


To derive physical amplitudes, we collect our results for the two- and three- point functions. 
For the renormalized two-point functions one finds the following expressions: 


('®) 


r 

^ XX 

B/ix 


B/ix 

Bh/i 





^ki 



+ 4(s) 



—y/Nv ^ 

(.+ 


1 _ 

^ N 

^B{s) 


N 


\ ^ 
) 


(32) 


The leading order contribution is contained in the momentum dependent mass. The 
GBs, labeled hy k,l E {1 ,... ,N — 1}, do not mix with the x, h fields. The wave function 
renormalization constants appear because of our definition of Bi. As mentioned above, 
we used only the wave function renormalization of the GBs, because we are interested in 
amplitudes without external Higgs particles. Thus it is not necessary to compute Z/j. With 

Qf _ J^/ 

the methods described in the previous section, it is easy to calculate Z^ — Zh = —. The 

inverse of the two point function matrix gives the propagators. As explained earlier, the 
diagonal GB propagators still have their pole at s = 0. One finds: 


Dhh{s) = 


with 


Z;, s - M\s) 

1 L 1 rm^(s) 


1 + 


M'^{s) m?{s) y N 

The other propagators are: 

Dhxjs) _ IZhDhh{s) 1 


A(s) -|- 2C{s) -\- 




B{s) 


-Pxx(g) 

The three-point functions are: 


Z^Dl,{s) Nl 
ZhDhh{s) 1 rm2(s) 


A(s) -|- C*(s) 




r/i7r7r(s) — 


2E{s) 


N 


- XTT’r 


(s) = -1- 


F{s) 

N 


(33) 


(34) 


(35) 


With these ingredients, one can express the amplitudes of the reactions // —> f f and 
// —> zz. Note that for any Hff vertex one gets a multiplicative factor ^ZhjZ^. The 
decomposition of the amplitudes into one-particle irreducible contributions is depicted in 
figure 
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Figure 8: Decomposition of the NLO 1/N corrections for the amplitude M{ff 
irreducible vertex eorrections and dressed propagators. 


zz) into 


Miff^f'D 

Miff^zz) 


1 


s-M^{s) 

1 m?{s) 

s - M‘^{s) y/Nv 


(36) 




The amplitudes can be measured directly at a future muon collider, and also can be 
related to Higgs production by gluon fusion at the LHC. The processes are available also to 
NNLO in usual perturbation theory, and thus are good quantities to compare the different 
theoretical approaches. One has simply to put = 4 and v = 123 GeV in our formulae. 

In figure ^ we plot the amplitudes of the two processes for several values of the coupling 
parameterized by mh, defined in eq. 0 - The Yukawa coupling is set to one. Note how 
the Higgs signal evolves from a sharp resonance to the damped regime, as the coupling is 
getting strong. 


We compare our results with the perturbative ones in figure 10 


The peaks of the NLO 1/N and NNLO perturbation theory result are aligned to allow 
for a physically meaningful comparison between the two calculations. For on-shell Higgs 
masses below 1000 GeV the curves practically coincide. From this figure it is obvious that 
the results of both approaches agree with each other up to the highest accessible energies of 
the LHC. Away from the resonance, the perturbative curves which we give are not expected 
to be a good approximation. Close to the vector boson pair production threshold, the 
equivalence theorem is not valid anymore, and the curves will be modified by contributions 
of order of the vector boson masses. 


As discussed in earlier publications |n.|, one finds in both approaches a saturation of 
the peak position. This is confirmed now also for the process // ^ zz. The saturation 
effect is shown in figure 

It turns out that the NNLO perturbative result allows a maximal peak position of about 
1100 GeV, whereas for the 1/N result one finds 975 GeV. In the case of // —> /'/' one 
finds 980 GeV for the perturbative and 930 GeV for the 1/N result. In this sense, the 
nonperturbative prediction is phenomenologically more restrictive. 
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5 Conclusions 


We described in some detail the computation of Higgs production processes in the nonper- 
turbative \/N expansion for processes which are relevant for future colliders. We want to 
stress that our calculation reproduces the perturbative results in the region up to around 
ruh = 800 GeV in a very precise manner, which shows the consistency between pertur¬ 
bative methods and the 1/N expansion including our tachyon subtraction scheme. Once 
the calculational tools are developed, the calculation is not more difficult than a two-loop 
calculation in the Higgs sector. Beyond rrih ~ 900 GeV, there are differences between the 
1/N expansion and perturbation theory concerning the position of the peak maxima, but 
the implications for Higgs-discovery physics are not too large up to ~ 1.1 TeV. This leads to 
the conclusion that the Higgs sector is theoretically well understood for the entire kinematic 
region relevant for Higgs searches at LHC [^. 

A saturation of the peak position at strong coupling seems to be a general feature of a 
heavy Higgs boson. 

Note that the 1/N approach disentangles the uncertainties from the renormalization 
scheme which is simply not present, from the cutoff effects due to the tachyonic pole. The 
effect of the tachyonic regularization becomes more important as the coupling is increased, 
modeling nondecoupling effects from a more fundamental theory containing the Higgs sector 
as an effective theory. 


Appendix 


The following formulae | 
graphs of eq. With its help the appearing two- 


IH] are used to get the integral representations of the 1/N 

and three-loop topologies can be cast 


into the form of one-loop integrals. 


Three-point massless scalar integral: 


with 




Li2(-—)+Li2(-—) 

^ Vi'' V 


+ 


87r2 FVA 


-Ul' 




'V2- 


ui ,2 = {l + b-a±VA)/2 
vi,2 = (1 — 6 -I- a ± a/A)/ 2 
A = 1 — 2(a + b) + {a — 6)^ 
a = p^/k^ 
b = q^/k^ 


(38) 
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Four-point massless scalar integral: 

We use the notation kij = the Mandelstam variables which are defined by 

the external momenta ,,, 4 ; 


Df){ki, k2, /cs) = Do{ki2, kis, k^, ^ 23 , A:24, A;34) = 

2 


1 


1 


-7;^os^i-xk) 


(47r)2a(xi - X 2 ) I 2 

k34 — i5 


-Li2(l + 
-Li2(l + 


/ci 3 - i5 
^24 — iS 


Xk) - r][-Xk, 




k34 — i6 
ki3 - iS 
^24 ~ 




ku-iS ) 'V kx2-i4 
+ log(-Xfe) [log(/ci 2 - iS) + log(/ci 3 - iS) 


ksi — i5 
ki3 - iS 
^24 ~ 

ki2 - iS 


Xk, 


Xk, 


- log(/ci4 - i5) - log(/C23 - i4) 


(39) 


whereby, 


a;i ,2 

a 

b 

c 

d 

V{u,v) 


—6 ± \/b‘^ — 4ac — Aiad5 
2a 

^24^34 

^13^24 + ki2k34 — kiik23 

ki2ki3 

k23 

log(uu) — log(tt) — log(u) 


The cut of the logarithm lies along the negative real axis. An infinitesimal replacement iS 
is included to stay on the correct Riemann sheet if an argument vanishes. 

Both formulae allow for analytical continuation of an external momentum which is 
important for our numerical purposes. This was checked by recalculating several known 
two-loop diagrams with the same methods as used for the computation of the 1 /N graphs. 
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